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QQ I Abstract 
(N 

■ We make a model-independent analysis of all available data that indicate neutrino 
i oscillations. Using probability diagrams, we confirm that a mass spectrum with two 

nearly degenerate pairs of neutrinos separated by a mass gap of ~ 1 eV is preferred over 

■ a spectrum with one mass eigenstate separated from the others. We derive some new 
relations among the four-neutrino mixing matrix elements. We design four- neutrino 

PLhI mass matrices with three active neutrinos and one sterile neutrino that naturally in- 

pi • corporate maximal oscillations of atmospheric and explain the solar neutrino and 

O , LSND results. The models allow either a large or small angle MSW or vacuum oscil- 

lation description of the solar neutrino deficit. The models predict (i) oscillations of 
either z^e — > Ur or Vs iu long-baseline experiments at L/E ^ 1 km/GeV, with 

^ . amplitude determined by the LSND oscillation amplitude and argument given by the 

^ \ atmospheric (5m^, and (ii) the equality of the Ue disappearance probability, the disap- 

pearance probability, and the LSND i/^ v^. appearance probability in short-baseline 
experiments. 



1 Introduction 



The long-standing solar neutrino deficit |T|, 0, the atmospheric neutrino anomaly P, H, H, H, 
and the results from the LSND experiment on v^, neutrinos from /i+ decay and v^. neutrinos 
from vr"*" decay can each be understood in terms of oscillations between two neutrino 
species |p. Interestingly, the solar, atmospheric, and terrestrial (LSND) neutrino oscillations 
have different L/E and therefore require different neutrino mass-squared differences 5im? 
to properly describe all features of the data. For example, if the atmospheric and LSND 
Sm? scales are the same one forfeits the recently reported zenith-angle dependence and 
up/down asymmetry of the atmospheric neutrino flux ^. Alternatively, if the solar 
and atmospheric 5m^ scales ||T0|] are the same, the reduction in the solar neutrino flux is 
energy-independent, contrary to the three solar experiments which infer different oscillation 
probabilities in different neutrino energy regions |]TT1. Since three distinct mass-squared 



differences cannot be constructed from just three neutrino masses, the collective data thus 
argue provocatively for more than three oscillating flavors. An alternative but less compelling 
possibility is to introduce new lepton-fiavor changing operators with coefficients small enough 



to evade present exclusion limits, but large enough to explain the small LSND amplitude [|T2 

If all of the existing observations are confirmed, a viable solution is to invoke one or more 
additional species of sterile light neutrino [jl3l, thereby introducing another independent 



mass scale to the theory. The additional neutrino must be sterile, i.e. without Standard 



Model gauge interactions, to be consistent with LEP measurements oi Z ^ vu |jTj]. The 



introduction of a sterile neutrino to complement the three active neutrinos has had some 



phenomenological success |[T5[| . 



In this paper we propose and study mass matrices for four-neutrino models (three ac- 
tive plus one sterile) that can accommodate all the present data. Once a fourth neutrino is 
admitted to the spectrum, it is no longer mandatory that the mix with the Vr at the atmo- 
spheric scale. The z/^ may instead mix with the sterile or with some linear combination 
of Vg and v^-. Similarly, the v^. may mix with a linear combination of z/^ and v^. 

At first sight the mixing of a sterile neutrino with active flavor neutrinos seems to be 
stringently constrained by Big Bang nucleosynthesis (BBN) physics. The bound 



5m^sin^2^ < 10"' eV^ (1) 
on the mass-squared difference 5m? and the mixing angle of the sterile neutrino was inferred 



to avoid thermal overpopulation of the "extra" sterile neutrino speciesflg]. However, there 
are significant caveats to this bound. One is the fact that some recent estimates of using 
higher abundances of ^He yield considerably weaker bounds [0. Another is that a small 
asymmetry [ui, — nv)/n^ ^ 7 x 10~^ of flavor neutrinos (but large compared to the present 
baryon asymmetry AnB/n^ ~ 10"^°) at t > 0.1 s is enough to suppress — Vg oscillations 



and then the bound of Eq. (^ does not apply |T8[. Such asymmetries, in fact, can be 



generated with the kind of model parameters considered herein (as shown in Ref. |19, |20| 
In light of this observation that BBN may allow sizeable mixing between sterile and active 
neutrinos, we consider both the small and large mixing with sterile neutrinos in this work. 

We review all existing data that indicate neutrino oscillations, and then perform a model- 
independent analysis of the data using four-neutrino unitarity constraints. A very useful 
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tool for this unitarity analysis is the set of probability rectangles, which we explain and 
exploit. We draw several model-independent conclusions for the four-neutrino universe. 

We design a five-parameter neutrino mass matrix which can account for each of the 
three viable solar solutions and accommodate the atmospheric and LSND observations. The 
three solar possibilities are the small-angle matter-enhanced (SAM) pi], E2], |23|, large- 



angle matter enhanced (LAM) |£5[ and large-angle vacuum long-wavelength (VLW) |£6[ 27 1 
explanations of the solar neutrino deficit. Our mass matrix yields maximal oscillations 
of atmospheric z/^. We consider the possibility that the solar data is explained by 
Ug or i/g — > z^T oscillations, in which case the atmospheric neutrino data is explained by 
either — > or z/^ — Us oscillations, respectively. We also consider the possibility that 
both atmospheric and solar neutrino oscillations have Ug and components. Lack of z/^-z/r 
discrimination in the present data is the major source of ambiguity in the four-neutrino 
model. We discuss how future experiments can resolve this ambiguity. 

In Sec. 2 we summarize the oscillation probability formulas and utilize a probability for- 
malism, based on unitarity of the mixing matrix, which permits a simple visual represention 
of mixing. In Sec. 3 we begin with a brief discussion of the three classes of experiments and 
the neutrino mass and mixing parameters needed to explain them. We then use probability 
rectangles to display the inferences from the data for any four-neutrino scheme. In Sec. 4 
we employ the probability rectangles to argue against a neutrino mass spectrum with one 
eigenstate separated from three other nearly-degenerate states (which we will refer to as the 
1+3 spectrum) in favor of two nearly degenerate mass pairs (which we will refer to as the 
2+2 spectrum). We also derive some new relations among elements of the mixing matrix 
that result from data and unitarity which are satisfied in a four-neutrino model for certain 
ranges of the parameters. Then in Sec. 5 we present a mass matrix whose eigenvalues consist 
of a nearly degenerate neutrino pair at ~ 1.4 eV and another nearly degenerate pair at low 
mass, as illustrated in Fig. 1. We show how the existing data almost uniquely fixes the model 
parameters (once a solar scenario is specified) and strictly determines what new phenomenol- 
ogy the model predicts. In Sec. 6 we derive expressions for the oscillation probabilities in 
our models in terms of the current neutrino experimental observables. We present the model 
predictions in Sec. 7. The new observable signature for the model is ^ Vt or z/e ^ z/^ 
oscillations for L/E ^ 1 km/GeV, depending on whether the atmospheric oscillations are 
Us, respectively. Section 8 contains some discussion, and a summary. 



2 Formalism 

2.1 Oscillation amplitudes 

To simplify the analysis of the available data, we will ignore possible CP violation and 
work with a real-valued mixing matrix U. Accordingly, the general formula for the vacuum 



oscillation probabilities becomes 28 



P{Ua Vp) = - 4 ^ UakUpkUajUpi sin^ Aj-fc , (2) 

k<j 

where Aju = 6m]^L/4E = 1.27{6m]JeV^){L/km)/{E/GeY), 5m]j^ = m] - ml and the 
sum is over all j and k subject to k < j. 
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For oscillations of two neutrinos, the oscillation amplitude (i.e., the coefficient of the 
sin^Ajfc term) is given by sin^2^, where 6 is the mixing angle between the two neutrino 
states. More generally for an arbitrary number of neutrinos, the amplitude of the Va to 
oscillation in the absence of CP violation is seen to be 

A'^^ = -AY,U^,UpjU^kUpi,, a ^(3, (3) 

where the sum is over mass states with mass-squared differences appropriate for the L/E 
of the particular experiment. We note that the oscillation amplitudes defined here are 
only for those oscillations at a particular A scale in Eq. (^). We will use subscript labels 
on the amplitude to identify the A scale (which is determined by the relevant and 
L/E) for the particular experiment: "sbl" will denote short-baseline experiments such as 
LSND, "atm" will denote atmospheric and long-baseline experiments, and "sun" will denote 
extraterrestrial experiments, especially those with solar neutrinos. 

We will use superscripts on the amplitude to identify the oscillation flavors, unless it is 
obvious from the context; in the absence of CP violation, = A^". With four neutrino 
states, [/ is a 4 X 4 mixing matrix. We also define the amplitude for Ua disappearance 

^a^ = ^ , (4) 

where /y represents a sum over neutrino flavor eigenstates other than Ua- The mixing- 
matrix elements U, and therefore the amplitudes A, depend on the environment, e.g., matter 
vs. vacuum. Throughout this paper we will quote values for the oscillation amplitudes in 
vacuum. 



2.2 Probability rectangles and a theorem 



The "probability rectangles" used by Liu and Smirnov visually illustrate the mixing of 
the flavor eigenstates among the mass eigenstates. To construct the probability rectangles, 
we introduce the notation 

Paj = \Uaj\ ; (5) 

such that Paj is the probability that the a*'^ flavor state is found in the j^^ mass state, or, 
alternatively, the probability that the j^^ mass state is contained in the a^^ flavor state. 
Therefore, when CP-violation is neglected, the real mixing-matrix elements are determined 



by the probabilities up to a sign: Uaj = ±^^JPaj- In principle, these signs may be determined 
by arranging for orthogonality of the rows, and columns, in the unitary mixing matrix U. 
By unitarity of U we have 

E^-. = l (6) 

a 

for each mass state j, and 

j 

for each flavor state a. Thus, if each mass state is represented as a rectangle of unit area, 
then the fractional area assigned to Paj within the rectangle is a graphical representation of 
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the value of Paj- The probabihties Paj depend on whether the environment is vacuum or 
matter. For consistency, we will always display vacuum probabilities in the rectangles. When 
the probability rectangles are displayed along a vertical axis labeled with mass-squared, the 
(5m^ values relevant for the various experiments are readily visualized. Figure 2 gives an 
example of the probability rectangles for a four-neutrino model. An inverted 2+2 mass 
spectrum, where the solar z/g oscillation is driven by the separation of the heavier two states 
and the atmospheric z/^ oscillation is driven by the separation of the lighter two states, may 
also be possible, but is not considered. 

The following mini-theorem will prove to be useful: 
In the absence of matter effects, the amplitude A°"^ is independent of how the P^j probabilities 
are partitioned among the mass eigenstates. 

The proof of this statement relies on the insertion of J^p^a UpjUpk = —UajUak into 

= -4E/3^aEi,fc>, U^jUakUpjU(3k, tO get 

A-^ = 4 5]P„,P„fc, (8) 

j>k 

where, as in Eq. (^), the sum in Eq. is over all mass states with mass-squared differences 
appropriate for the L/E of the particular experiment. In Eq. (H), Aagt is manifestly indepen- 
dent of the partitioning of the P^j probabilities since it involves only the Paj- This theorem 
demonstrates the limitations on information derivable from disappearance experiments. 

The minitheorem fails in the presence of matter effects because the partitioning of the 
flavor probabilities including Paj are altered. That is, with matter effects the amplitude Aa^^i 
does depend on how the P^j are partitioned among the mass states. Matter will also alter 
the oscillation wavelength, causing further changes in the phenomenology of experiments 
sensitive to the oscillations rather than their averages. Matter effects have the potential to 
resolve the Us-Ur ambiguity, as do some other measurements. We discuss these possibilities 
in Sec. 7. 



3 Experimental constraints 

3.1 Short baseline: LSND, reactors, and accelerators 

The LSND experiment reports positive appearance results for z/^ — > z/g oscillations from 
fi^ decay at rest (DAR) and for z/^ z/g oscillations from vr"*" decay in flight (DIE). The 
DAR data has higher statistics, but the allowed regions for the two processes are in good 
agreement. There are also restrictions from the null results of the BNL 13-776 [^] and 



z/p 



KARMEN z/^ z/g oscillation search experiments. The combined data suggest 
vacuum oscillation parameters that lie approximately along the line segment described by 

0.3 eV^ < Sml, = ^^^^ < 2.0 eV^ . (9) 

y^sbi) 

However, values for Sm^f^i as high as 10 eV^ are also allowed for A'^^i ~ .0025, although values 
above 3 eV^ are disfavored by the r-process mechanism of heavy element nucleosynthesis in 



supernovae |3l|. 
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There are also relevant data from the Bugey reactor experiment which searches for z/g 
disappearance and from the CDHS|Q and CCFR[0 experiments which set bounds 



on disappearance. 

The combined short baseline data set for A'^^, A^^'^, and A'^^ will be used in Sec. 4 to argue 
against a hierarchical neutrino mass spectrum in favor of two pairs of nearly degenerate 
masses in the four-neutrino spectrum. 

3.2 Atmospheric data 

The atmospheric neutrino experiments measure z/^ and (and their antineutrinos) created 
when cosmic rays interact with the Earth's atmosphere. One expects roughly twice as 
many muon neutrinos as electron neutrinos from the resulting cascade of pion and other 
meson decays. Several experiments || obtain a z/^/z/g ratio that is about 0.6 of the value 
expected from detailed theoretical calculations of the flux ||3^. The Super-Kamiokande 



(SuperK) experiment has collected the most data and analysis [H indicates that their results 
for contained events can be explained as z/^ — > Ur oscillations with 



3 X 10^^ eV^ < 6ml^^ < 7 x 10"^ eV^ , 0.8 < A'^t, < 1.0 . (10) 

The high end of each range is favored. 

Independent of flux normalization considerations, the — > z/g oscillation channel is 
strongly disfavored by the zenith angle distributions of the data and by the up/down 
asymmetry separated into "muon-like" (z^^) and "electron-like" (z/g) events 0, which yield 
an up-to-down ratio of 0.52lg:[]^ ± 0.01 for /x-hke events and 0.84l|j;}^ ± 0.02 for e-like 
events (the expected values are close to unity). Furthermore, the recent CHOOZ z/g dis- 
appearance experiment excludes z/g — > z/^ oscillations with large mixing A^f^ ^ 0.2 for 
Sml^ > lO-^eV^ 0. 



In a four-neutrino context, another possibility for the atmospheric neutrino oscillations 
is z/^ — > Us- Oscillations of this type in principle could be affected by matter due to the 
different neutral current interactions of and Ug. However, for the contained events (with 
lower energy) these effects are small, especially for larger values of hence, the 



allowed regions for z/^ Us should be similar to those for z/^ — >■ z/,-. For events at higher 
energies the matter effects could begin to be appreciable; a definitive test requires more data. 

3.3 Solar data 

The solar neutrino experiments 0] measure z/g created in the sun. There are three types 
of experiments, z/g capture in CI in the Homestake mine, z/g — e scattering at Kamiokande 
and Super-Kamiokande, and z/g capture in Ga at SAGE and GALLEX; each is sensitive 
to different ranges of the solar neutrino spectrum and measures a suppression from the 
expectations of the standard solar model (SSM)[|l|. 

For z/g — >■ Us oscillations in the sun (in which case atmospheric neutrino oscillations are 
1^1,. i^T in our model) the allowed parameter ranges at 95% C.L. for the small-angle 



matter-enhanced solution are given in Table ^ The solution is based on the SSM fluxes 
in Ref. [^. Approximate parameters for the large-angle matter-enhanced |^ and vacuum 
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long- wavelength solutions [J0| for z/g oscillations of solar neutrinos are also shown in 

Table 0. If the solar neutrino deficit is caused instead by Ue oscillations (and the 

atmospheric oscillations are — > Ug), then the allowed solar parameter ranges for the three 
solar cases are slightly different ; see Table |l|. The exact values of the parameters may 
change as new data from SuperK become available and when fits are made with the new 
solar flux calculations. 

In any of the matter-enhanced scenarios it is also necessary that the eigenmass mi 
associated predominantly with z/g be lighter than the eigenmass mo associated predominantly 
with the neutrino into which the z/g is oscillating (i.e., or z/^), so that it is rather than 
z/e that is resonant in the sun. For the vacuum solutions the ordering of mg and mi does not 
matter. Alternate scenarios where the z/g is predominantly associated with the heavier two 
states and z/^ is predominantly associated with the lighter two states are also viable. 

For z/g Vr oscillations in the two-neutrino approximation the propagation equation for 
the neutrino states in the charge-current basis is pll, H2|, E3| 



d 
I— 
dt 



ye 



}_ ( A^/2GFENe 5m^ sin 2d 
AE \ 6m^ sm2e 26m^ cos2e 



(11) 



where A'^g is the electron number density. For z/g 
is instead 



Vg oscillations the propagation equation 




}_ ( 4:V2GfE{N, - ^Nn) Sm'^ sin 29 
4E \ Sm"^ sm2e 2Sm'^cos29 




(12) 



where Nn is the neutron number density. For the small-angle matter-enhanced case the non- 
adiabatic approximate solution for neutrino propagation is appropriate and the oscillation 
probability for a neutrino of energy E^, is 



sun J 



where in this case ^ labels either r or s, and 



exp 



4:E„ 



Asun {d\ogN/dL), 



(13) 



(14) 



is the Landau-Zener transition probability and is either Ne (for z/g — > z/^ oscillations) 
or A'g — ^Nn (for z/g Ug oscillations). The quantity {d\ogN/dL)c is the appropriate 

logarithmic density gradient in the sun at A^^^** = 5im?g^^\J 1 — Altn/ {2\/2GFEy)^ the critical 
density where maximal oscillations (resonance) occur. For the large-angle matter-enhanced 
case, the neutrino propagation is adiabatic and 



(15) 



assuming the neutrinos are created where the electron density is well above the critical 
density. For the vacuum long-wavelength solution the oscillation probability is just given by 
the usual vacuum expressions. 
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3.4 Oscillation lengths and amplitudes summarized 



In neutrino oscillation descriptions of the solar, atmospheric, and LSND data, a distinct 
oscillation wavelength and oscillation amplitude is required for each of the three data sets. 
Experimental uncertainties allow for some latitude in these amplitudes and wavelengths, 
and for the solar data, there are three isolated islands of viability in the (5m^-amplitude 
plane see Fig. 3. The day-night asymmetry measurement, found to be small in the 



recent SuperK data, removed about half of the previously viable solar regions |45 



The vacuum oscillation wavelength is linear in the neutrino energy, allowing further pos- 
sibilities that are summarized in Table The chosen neutrino energies are typical for solar, 
and reactor sources (5 MeV), pion facilities (100 MeV), for contained (2 GeV), partially- 
contained (10 GeV), and throughgoing (100 GeV) neutrino events in underground detectors, 
and for astrophysical sources (1 TeV). Although full oscillation wavelengths are also listed in 
Table ^ oscillation effects may well be measurable for a fraction of an oscillation wavelength 
or as an average over many oscillation wavelengths. Throughgoing and partially contained 
atmospheric neutrinos may show nodes as a function of L/E. Further possibilities arise when 
the matter effect of the earth is included in the oscillation physics. We consider earth-matter 
effects in Sec. 7.6. 



3.5 Inferences from data 

We consider first the probability rectangles for the atmospheric and CHOOZ data. The 
atmospheric data indicate Sm^^^ ~ 5 x 10~^ and nearly maximal flavor-changing mixing of 
with u^. The present data do not distinguish between u^- or Ug as the dominant state into 
which Ufj, mixes. The probability rectangles for the atmospheric scale are displayed in Fig. 4a. 
We label the two masses defining the atmospheric scale as 1^2 and z/3, with = 5^32. 

Because of the ambiguity we show the union P^+Ps rather than the partitions into Pr 

and Pg. For maximal z/^-z/^ mixing, one must choose P^2 ~ -P/^s ~ Pt2+Ps2 ~ Pr3+Ps3 ~ 1/2. 

Next we consider the pair of mass eigenstates whose mass-squared difference is fixed by 
the solar scale. We provisionally investigate a four-neutrino mass spectrum that consists 
of two pairs of nearly degenerate neutrinos separated by the LSND scale ^m^^; ~ eV^ (to 
explain the LSND result in terms of oscillations). We argue in Sec. 4 that the data favor this 
spectrum over a spectrum with one mass separated from three relatively degenerate masses. 
We label the second pair of mass states as z/q and z/i, and define Smj.^^ = 5m^g. Since P^2 
and Pfj,3 sum to near unity, P^o and P^i must be small. Thus the probability rectangles 
for the z/q and z/i states appear as shown in Fig. 4a. Accordingly, the LSND amplitude for 

l^f, l^e, 

^sM = -4 [U,3Ue3 + U,2Ue2] [U.^U,, + U^oUeo] , (16) 

is necessarily small. We emphasize that the smallness of the LSND z/g-z/^ mixing is an 
inevitable consequence of the large mixing of z/^ to at the atmospheric scale and the con- 
straints of unitarity, independent of particular model considerations including rearrangement 
of the neutrino mass spectrum. 

Four-neutrino unitarity may be used to rewrite eqn. ( [TE| ) as 

^sbl = 4 K3Ue3 + U^2Ue2? = 4 \U,^U,^ + U^^UeO?- (17) 
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Written this way, it is clear that the LSND data is bhnd to the partitioning of u-r and Vs 
in the probabihty rectangles of mass states z/q and Vi. This flavor ambiguity is shown in 
Fig. 4a. 

With the identification = 5m^Q, we may use Eq. (P) to write the solar z/g- 

disappearance amplitude as 

A± = AP,oPei. (18) 

Because matter in the sun may exert a significant effect on propagating neutrinos, the values 
of Peo and Pel for the sun have some sensitivity to the state v^. or z/^ into which z/g oscillates. 
However, the sensitivity of present data to this difference is weak, and there is considerable 
freedom in assigning or z/^ or a linear combination thereof as the mixing partner to Ve- This 
z/^-z/g ambiguity for Ve mixing at the solar scale is complementary to the Vg-Vr ambiguity for 
z/^ mixing at the atmospheric scale. Potential measurements to resolve the z/^-z/t- ambiguity 
at the solar scale will be discussed in Sec. 7. 



3.6 The three solar solutions 

The Peo and Pgi partitioning specifies whether the solar model is a small-angle model or 
a large-angle model. As can be inferred from Eq. (0), with nearly-equal partitioning of 
Peo and Pel, the mixing amplitude is near maximal (large angle). With highly nonequal 
partitioning, i.e., Peo -C Pei or Pei -C Peo, the mixing amplitude is small. Of the three 
viable solar neutrino options, SAM falls into the small angle category, while LAM and VLW 
fall into the large angle category. The probability rectangles for the small and large angle 
classes of models are shown in Figs. 4b and 4c. Recall that in order to obtain the MSW 
resonant enhancement required for the SAM and LAM solutions, it is necessary that the 
state which is predominantly Ve be the lighter of the two mass states, z/i. Qualitatively, 
LAM and VLW are distinguishable in their probability rectangles only by the choice of value 
for 5m1^^. Quantitatively, the two solutions and the VLW solution are distinguishable in 
ways which are discussed in Sec. 7. 

If the active-sterile mixing is small, then all ambiguities in the probability rectangles 
are resolved: the large atmospheric mixing must be y^-Vri and the solar solution must be 
small-angle SAM with Ve-Vs mixing. The probability rectangles for this model are shown 
in Fig. 2. This particular solution has recently been analyzed in the context of the minimal 
four-neutrino mass matrix EBI, B7 . 



4 Mass spectra 

4.1 Argument against a 1+3 mass spectrum 

It has been shown by Bilenky, Giunti and Grimus that a hierarchical ordering of the 
four-neutrino spectrum (implying one dominant mass) is disfavored by the data when the 
null results of reactor and accelerator disappearance experiments are included. We will refer 
to this spectrum as the 1+3 spectrum, defined as one heavier mass state separated from 
three lighter, nearly-degenerate states, or vice versa. We demonstrate the argument with a 
set of logical steps similar to theirs. 
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Assume a mass spectrum with one heavy mass well separated from three other nearly- 
degenerate states and let the heavy mass state be labeled as u^. Then the LSND mass-squared 
scale is Smlf^i ~ — '^^31 — ^^30 LSND amplitude is 

A':',, = -4 U^sU^s [E U,jU,,] = 4 P,,P^s, with P,3 + P>.3<1. (19) 

On the other hand, the z/g and z/^ disappearance experiments at reactors and accelerators are 
also sensitive to the LSND scale. These experiments measure the disappearance amplitudes 

A% = 4 Pes E Pej] = 4 Pes [1 " Pes] (20) 

and 

A':^, = AP,,[1-P,,]. (21) 

The second equalities in Eqs. ([I9| ) and (|20| ) (see Eq. (|^)) follow from unitarity of the mixing 
matrix. The three amplitudes in Eqs. (|T9| ) -(l2Tl ) depend on just two parameters, and so are 
interrelated. All three of these amplitudes are constrained by experiments to be small. A 
priori then, Pe3 and P^s may both be small, or one (but not both) may be near unity with 
the other small. The fact that A'^^i is an appearance observation rather than a bound means 
that if Pes and P^3 are both small, they cannot be too small. 

In the 1-1-3 model, the atmospheric scale does not involve the heavy state u^. Without 
loss of generality we label the state which determines the atmospheric scale as i>2. Then from 
Eq. (H) the atmospheric z/^ disappearance oscillation amplitude is given by 

A>:L = 4P^2(Pmo + P^.l) < (1 - Pf.sf , (22) 

where the inequality comes from maximizing 4:P^2{PiJ.o + Pfii) subject to the constraint P^o + 
P^i-|-P^2 = 1 — -P^iS- The SuperK data indicate that z/^ is maximally mixed at the Sni^^^^ scale, 
i.e., there is little z/^-content available to the z/3 state. Quantitatively we have A^^f^^ ^ O-S^ 
which implies P^a < 0.11. Since P^a is small, Eq. ( PT| ) becomes 



A^,^^4P^3«1. (23) 

The probability rectangles for the 1-1-3 model with small P^s are presented in Fig. 5. Note 
that it is the zenith-angle, or up/down asymmetry data, which really establishes Sml^^^ as 



different from Snr^f^i, that is crucial for the argument |^ 



We are left with the possibilities of Pes being small or near unity. As can be seen in Fig. 5, 
if Pe3 is near unity, then there is little Pe^ to distribute over the three lighter mass states. 
In particular, the solar amplitude A^^^ = iP^oPei, where the solar scale is 5m1^^ = StuIq, is 
second order in small quantities, too small for even the SAM solution {Asam > 2.5 x 10"^) 
to the solar flux. This may be easily quantified. If Pes were near unity, we would have 

A%c^4[l-Pes]«l. (24) 

Together with unitarity, this in turn bounds the magnitude of the solar amplitude: 

Attn = ^PeoPel < (1 " Pes)' ^ ^{^1)' , (25) 



10 



where the inequahty in Eq. (|25|) comes from maximizing 4Peo-Pei subject to the constraint 
Peo + Pel < 1 — Pe3- The experimental upper hmit on A^fi from the BUGEY experiment ^2 



is about 0.1 for 5m^j; < 2 eV^, which disallows even the small-angle solar solution. We 
conclude that Pes is small, in which case 

A:ti^AP,,«l. (26) 



Thus, both Pe3 and P^s must be small in the 1+3 model, and from Eqs. ([19|), (P3|), and 
^), we infer the relation 

^sbl — ^^sbl ^sbi- l^' J 



However, the experimental upper bounds on the disappearance amplitudes A^li and 
^sbi IH and the measured appearance result for A'^^^ M are not compatible with Eq. (|2 



sbl ^'-'-'^ '^^^^ liiCdDuicLi appeal icDuit lui 

thereby disfavoring the 1+3 model. For example, for Sm'^.f^i = 0.3 eV^, A^fi^ < 0.035 from 
Bugey, y4^^ < 0.8 from CDHS, which implies A'^^^ < 0.007; however, for this value of Sml^i, 
the LSND data indicate A^^^i > 0.04. The LSND results are presented in terms of maximum 
likelihood rather than confidence level limits, so it is not straightforward to state an exclusion 
probability. 

Put another way, A^^^ is large enough that the Bugey and CDHS limits force one of 
Pe3 and P^3 to be small and the other to be large, but this is ruled out by the solar and 
atmospheric data. The constraints on Pes and P^s from the three short-baseline amplitudes 
^Ibh ^Ibh ^sbh ^^e atmospheric amplitude A^tm and the solar amplitude A^^^ (from 
Eqs. ([T9|) , (pO|) , (^Tl) , (p2|), and (|25|) , respectively) are conveniently summarized in Fig. 6 for 



two different values of 5m^j,;. 

The measured values and bounds for the short-baseline appearance and disappearance 
amplitudes depend on the magnitude of 5m ^j,;. (There is effectively a suppressed 5m ^j,; third 
axis in our Fig. 6, which samples only two particular values of 5m ^^;.) For certain allowed 
values of 5m^ (e.g., at 1.7 eV^ and 0.25 eV^ according to Fig. 2 of |^) the violation of 
Eq. (pTl) is mild, and the 1+3 model is just barely incompatible with the data; see, e.g.. 
Fig. 6a. 

The argument against the 1+3 model does not depend on the sign of 5m.^^;. This means 
that the inverted 3+1 model with the three nearly degenerate mass states heavier than the 
remaining state is equally disfavored. 



4.2 2+2 mass spectrum 

We now turn to the favored class of four-neutrino models, namely those with two nearly 
degenerate mass pairs separated by the LSND scale as displayed in Fig. 1. It is interesting to 
see how this "pair of pairs" mass spectrum of four-neutrino models realizes the dependency 
among A'^^i, A^^, and A^f,^ which conflicted with the 1+3 model. Let and z/i label the pair 
of the nearly-degenerate mass eigenstates responsible for the solar oscillations, and z/2 and 
Us label the pair of the nearly-degenerate mass eigenstates responsible for the atmospheric 
oscillations. 

The expressions for the oscillation amplitudes are 

A>:',i = 4 \U,2U^2 + U,sU,3\' = 4 \U,oU,o + U.^U.^W (28) 
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= 4(Pe3^el + Pe^PeO + Pe2Pel + PelPeo) = 4 (1 - a,), (29) 

and 

A^^ti = 4(P^3Pmi + ^m3^/.o + Pm2^a«i + ^/.2Pmo) = 4 (1 - (t^) , (30) 

with 

<ya = \Ua2? +\Ua3? = Pa2 + Pa3. (31) 

The Schwartz vector inequahty \ve ■ v^\'^ < \ve\'^ \v^\'^ apphed to the vectors = {Ue2,Ue3) 
and = {Uf,2, U^z) then gives 

A:i,<Aa,a,. (32) 
Furthermore, in the 2+2 model the solar oscillation amplitude is 

= 4PeoPel < (1 - CTe)' , (33) 

and the atmospheric oscillation amplitude is 

^± = 4P^2Pm3<^J, (34) 



where the inequalities in Eqs. (0) and (|3^) come from maximizing the expressions subject 
to the constraints Peo + Pel = 1 — cTe and P^2 + P^z = respectively. 

If the vector inequahty in Eq. (^2]) is saturated, then A^J^i, Alfi, and each has the 
same functional dependence on two parameters as it did in the 1+3 model (tXe has replaced 
Pe3 and cr^ has replaced P^s). Then the previous argument that the LSND, Bugey and 
CDHS data require one parameter to be small (<^ 1) and the other large (~ 1) applies. The 
argument is unaffected if the vector inequality is not saturated. As before in the 1+3 case, 
the solar constraint indicates that Ue must be small. This time however, unlike the 1+3 
case, the atmospheric constraint involves cx^ and not (1 — o"^)^, and can be met if is large 
(~ 1). Therefore the constraints of the data can be satisfied by assigning Ue dominantly to 
one pair of mass states and z/^ dominantly to the other pair. Instead of Eq. (^) pertinent 
to the 1+3 spectrum, we obtain for the 2+2 spectrum 

ATu < Kti ■ (35) 

This bound is linear in the small disappearance amplitudes, and is easily satisfied by the 
data. For example, the tightest constraint on A^^^ is about 0.02 for (5m^^; = 0.6 eV^, while 
the LSND data indicate A^li can be as low as 0.009 for this value of Sm^f^i. 

In Fig. 7 we have drawn the o"e-o"^ plot for the 2+2 model, analogous to the Pgs-P^s plot 
for the 1+3 model, for Sm^f^i = 1.7 eV^. The allowed regions with near unity (implying 
near-maximal mixing of z/^ in the z/2-t'3 pair) and cTg small (implying almost no mixing of z/g 
into the Z/2-Z/3 pair) show that the 2+2 model can comfortably accommodate the data. 

Since only mass-squared differences are important for oscillations, the inverted 2+2 
model, where the solar oscillations are driven by the mass-squared difference of the up- 
per mass pair and the atmospheric oscillations are driven by the mass-squared difference of 
the lower mass pair, is equally viable. 
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4.3 New results 



Two features of the data are especially noteworthy. The first is the remarkably high degree 
of isolation of Ve into one mass pair and into the other mass pair, as inferred from the 
bounds on the disappearance amplitudes. The second is the near saturation of the vector 
inequality in Eq. ( ^2]) by the LSND appearance amplitude A^^^ for 5m^^; ~ 0.3 eV^. 

Equations (^) and (pOf ) bound the degree to which z/e and z/^ are found in opposite 
pairs of mass eigenstates. Without loss of generality we assume that z/e is predominantly 
associated with z/q and z/i, and that z/^ is predominantly associated with z/2 and z/3. Then 
from the Bugey and CDHS data we find the constraints 

cTe ~ ^ A% < 0.016 (0.009) , (36) 

and 

l-^M^^^^6l< 0.013 (0.2), (37) 
respectively, for = 2 eV^ (0.3 eV^). We then deduce 

Ali < 0.065 (0.04) , (38) 
which can be compared to the LSND data 

A^fc' ~ 0.0025 (0.04) , (39) 



for these two values of ^m^^;. The near-saturation of the inequality in Eq. (|35| ) for 5m^^; = 
0.3 eV^ has very interesting implications. It means that is nearly parallel or antiparallel 

Ue2/U,^\ ~ |f/^2/t/^3| • (40) 



to -0^, which in turn indicates that 



This is a new result. 

Furthermore, the SuperK data suggest that z/^ is maximally mixed in the mass pair with 

atm 1 



mass-squared difference so for this pair, called 1^2 and z/3, that 



\U,2\ ^ \U,,\ ^ (41) 

Then Eq. (|40| ) implies 

\Ue2\ ^ \Uez\. (42) 

This is also a new result. In summary, if the oscillation parameters are indeed near the 
limits of the Bugey bound, the four-neutrino mixing matrix in the 2+2 model must satisfy 
Eq. (^Op, which implies Eq. (^) if the atmospheric z/^ mixing is maximal. 

We can derive additional constraints by considering a'^ = Peo + Pel and = P^o + -P/^i 
rather than ae and a^. The data requires a'^ to be large (~ 1) and cr^ small, and the Schwartz 
inequality reduces to 

ATm<A''J,, (43) 
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where the CDHS bound is 

A%^i < 0.05 (0.8) , (44) 

at Smlf^i = 2.0 eV^ (0.3 eV^). Because the inequahty in Eq. (^31) is not saturated by the 
data for any Srn^fji, a relation similar to Eq. ( ^0[ ) for Ueo, Uei, U^o, and U^i is not required 
in the 2+2 modeL However, the exphcit mass matrices we consider do have such additional 
relations; see Sec. 5. 

Finally, we mention a curiosity in the data which occurs for the pair of pairs mass 
spectrum with the matter-enhanced solar solutions (SAM and LAM). The linear mass split- 
ting at the heavier pair is rris — m2 ~ 5777,32/2777,3. If this pair is associated with the at- 
mospheric scale, we have ma — 7772 ~ 2.5 x 10~^eV (5m^j„/5 x 10~^eV^)(r?73/eV)~^. On 
the other hand, if the lighter mass pair is associated with the matter-enhanced solar scale, 
and ??7o ^ ??7i, then the linear mass-splitting of this pair is mo — mi ~ itiq ~ ^ ^^lun — 

2.5 X 10~^eV (5m^^^/10~^eV^)^. The two linear mass splittings within the pairs are nearly 
identical. While squared masses enter into the oscillation formulae for relativistic neutrinos, 
the more fundamental constructs of field theory, such as the Lagrangian and the resulting 
equations of motion, are linear in fermion masses (and quadratic in boson masses, these 
powers of mass being related to the dimensionality of the fermion field vs. the boson field). 
Thus it is a worthy enterprise to attempt to deduce linear neutrino-mass relations whenever 
possible. 



5 Mass matrix ansatzes 



5.1 Solar Up 



oscillations 



To describe the above oscillation phenomena in the scenario where the solar neutrino deficit 



is described by Ve oscillations and the atmospheric data by z/^ 

neutrino mass matrix ansatz 



we consider the 



M 



m 





£2 





\ 


















64 


1 


I 


es 


1 


£4 / 



(45) 



presented in the (i/^, z/g, z/^, i^r) basis (i.e. the basis where the charged lepton mass matrix 
is diagonal). By considering the field redefinitions \[' — \[' and \E' — > 75\E' one realizes 
that m, and at least one of ei and €4, and at least one of €2 and 63, may be taken as 
positive; we will take m and all four to be positive for simplicity. The mass matrix M 
contains five parameters (m, ei, 62, €3, €4), just enough to incorporate the required three 
mass-squared differences and the oscillation amplitudes for solar and LSND neutrinos. The 
large amplitude for atmospheric oscillations does not require a sixth parameter in our model 
because the structure of the mass matrix naturally gives maximal mixing of with (or 
with Vg if Vj. and Vg are interchanged). 

For simplicity, we have taken the mass matrix to be real and symmetric. The choice of a 
symmetric neutrino mass matrix is well-motivated in the context of oscillations, for what is 
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measured in neutrino oscillations are the differences of squared masses, which are eigenvalues 
of the hermitian matrix MM^ , which is itself symmetric when CP conservation is assumed. 
M is diagonalized by an orthogonal matrix U (real) and there is no CP violation. The 
are assumed to be small compared to unity, but not all necessarily of the same order of 
magnitude. The zero terms in the mass matrix could be taken as nonzero without changing 
the phenomenology discussed here as long as they are small compared to the terms shown. 
Also, the My^y^ term could be chosen different from 64 while still giving maximal mixing of 

and Vr since maximal mixing results from the large value of the M^^y^ matrix element 
relative to the diagonal M^^y^ and My^y^ elements, without any need for fine tuning of the 
difference \My^y^ — My^y^ \ . Here we choose to take the minimal form for M needed to describe 
the data and then derive the associated consequences. 

To a good approximation, the two large eigenvalues of the mass matrix in Eq. (^) are 

"^2,3 = Tm(lTtA + \tfj . (46) 

The values of the two small eigenvalues depend on the hierarchy of the ej. For the three 
solar cases we have: 



SAM 
LAM 
VLW 



62 < ei,e4 < es < 1 , (47) 
€1,62,64 < 63 < 1 , (48) 
61 < 62 < 64 < 63 < 1 . (49) 



The two small eigenvalues are then approximately given by 

SAM: mo ~ m6i , mi ~ m(6|64 — 6|/6i) , (50) 

(51) 



LAM : mo,i ~ f 



61 ± Jej + 46i 



VLW : mo,i - m [±62 + 6i64/2] . (52) 

These approximate expressions for the eigenvalues have been obtained by multiplying each 
6j by powers of a hypothetical parameter 6, where the number of powers of 6 assigned to 
each 6j depends upon the ordering in Eqs. (^)- (|49|) . For example, in the SAM case 63 are 
multiplied by 5, 61 and 64 by 5^, and 62 by 6^. Then each eigenvalue is written as an expansion 
in powers of 6, the coefficients of which may be solved for by requiring that the expression 
Ili{X — Xi) reproduces the eigenvalue equation for the mass matrix order by order in 6. Once 
the coefficients are found, S is set equal to unity. 

The eigenvalues in all cases have the desired hierarchy mi < mo -C 1712, m^, which 
gives the mass spectrum of the 2+2 model described in Sec. 4.2 and depicted in Fig. 1. 
The small relative mass splitting of the heavier masses m2,m3 is governed entirely by the 
parameter 64: Sml2 — 4m^64. The LSND —>■ oscillations are driven by the scale 
5m2i — Sm'^i ~ 5m2o — ^^lo — "^^j the atmospheric oscillations are determined by (Jmgj, 
and the solar z/g — > Ug oscillations are determined by SttiIq, the approximate expression for 
which can be obtained by Eqs. (^)-(|5^). The charged- current eigenstates are approximately 
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related to the mass eigenstates by 
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(53) 



where tan2^^ = 2e2/ei. The dots indicate nonzero terms that are much smaller than the 
terms shown. It is their smallness that suppresses mixing between and Vg. The mixing 
matrix V depends on just three of the original five parameters; it is independent of €4 and 
the overall mass-scale parameter m. Note that and couple predominantly to and 
z/g. The nearly-degenerate z/2 and 1/3 are seen to consist primarily of nearly equal mixtures 
of and v^. These results, illustrated in Fig. 2, conform to the qualitative arguments of 
Sec. 3 based on probability rectangles. 

It is noted that this mixing matrix not only satisfies the approximate equalities of 
Eqs. (^OD-(^), but in fact replaces the approximate equalities, derived from parameter- 
independent arguments, with exact equalities to first order in e^. Inspection of the mixing 
matrix reveals that our model predicts saturation of Eqs. (|35|) and (^31) to this order, i.e., 
A^^^i = A'^^i = A^li. A small improvement in the measurement of A^f^ or a modest improve- 
ment in the measurement of is predicted to show a positive disappearance signal. 



5.2 Solar Vp 



oscillations 



Another scenario, with solar z/g and atmospheric z/^ 

obtained by interchanging —>■ Ug and Ug 
basis is then 

/ £4 



Ug oscillations, is readily 
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m 
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-£2 



—Ut-. The mass matrix in the (z/^, z/g, z/^, z/,-) 

(54) 
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The eigenvalues and parameter hierarchies are still given by Eqs. (|46|)-(|5^). The mixing 
matrix is then given by 
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(55) 



where again tan 26' = 2e2/ei. 

In the VLW case, the parameter ei is negligibly small if the solar oscillations are maximal, 
and can be taken as zero without affecting the phenomenology. If this is done, then reference 
to the mass matrix shows that both z/g and derive their masses entirely from flavor non- 
diagonal couplings, and they are maximally mixed (analogous to the Vpri^s system). Also, 
if ei is taken as zero, then there are only four independent parameters needed in the mass 
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matrix, and just two in the mixing matrix. The derived 6 parameter becomes and the 
mixing matrix becomes very simple: 



U 



1 

V2 



( 



-1 M 

±1 1 63 63 

^63 -63 1 1 

V -1 ±1 / 



(56) 



5.3 Solar Up 



Vg and Ve 



Vr oscillations 



A more general scenario which is a mixture of the previous two is for solar neutrinos to 
undergo both z/g Vg and Ve oscillations. This is easily parameterized by replacing 

the Ug and z/^ states in Eq. with the rotated states z/^ and u'j. and defining 



cos a — sm a 
sin a cos a 



z/.« 



M 



m 



\ J^r / V sma cos a J \ i^- 

Then the mass matrix in the (z/^, z/g, z/^, z/^) basis becomes 

/ 61 cos^ a + 64 sin^ a 62 cos a + 63 sin a sin a (64 — 61) sin a cos a \ 

62 cos a + 63 sin a 63 cos a — 62 sin a 

sin a 64 cos a 

V (64 — 61) sin a cos a 63 cos a — 62 sin a cos a 64 cos^ a + 61 sin^ a / 

and the matrix which diagonalizes M is 

/ cos 6 cos a — sin 6 cos a 
sin 6 cos 6 

—63 sin 6* —63 cos 6* 



(57) 



(5^ 
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where tan 2^ = 262 /ei as before 
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(59) 



6 Oscillation probabilities 

6.1 Expressions for any baseline 

For the mixing in Eq. (^) (when the solar oscillations are z/g 



oscillation probabilities obtained from Eq. 
amplitudes smaller than 0{e^), are 



Us), the off-diagonal vacuum 
to leading order in 6,- for each Aj,- and ignoring 



P(iy„. 



2 cos^ 6'(sin^ A21 + sin^ A31) - sin^ A32 
+ 2 sin^ 6'(sin^ A20 + sin^ A30) - sin^ 29 sin^ Aqi 
63 sin^ A32 , 
sin^ 29 sin^ Aqi , 
sin^ A32 , 
sin^ 26*63 sin^ Aqi , 



(60) 

(61) 
(62) 
(63) 
(64) 



17 



where Aqi ^ A32 ^ A20 — A30 — A21 — A31 due to the neutrino mass spectrum. In our 
model, only the z/^-z/^ oscillation is suppressed beyond C(e^). 



6.2 Short baseline 

For small L/E only the leading oscillations A20 — A21 — A30 ~ A31 contribute, and the 
only appreciable oscillation probability is 

P{ye ^ y^,) ^ 4e^ sin^ A , (65) 

where A = rn^L/AE. From Eq. (^) we can fix two model parameters 

Sml, = m\ A^J, = A4. (66) 

Since the only short-baseline oscillation is z/g ^ z/^, these models predict the equality of 
the z/g disappearance probability, the z/^ disappearance probability, and the LSND z/^ —>■ 
appearance probability in short-baseline experiments. 



6.3 Long baseline 

For LjE typical to atmospheric or long baseline neutrino experiments, the oscillations in 
A assume their average values. The A32 oscillation is now evident, and the non-negligible 
oscillation probabilities in vacuum are 

P[v^ ^ Ur) ~ sin^ A32 , (67) 

P(z/, ^z/^)~e2(2-sin2A32) , (68) 

P(z/e ^ Ur) ~ sin^ A32 . (69) 

From Eq. (|6^) 

^^Itm = ^"^32 - 4m2e4 , Ai^i, = 1 , (70) 

which determines another parameter of the model. The model automatically gives maximal 
oscillations for atmospheric z/^'s, while oscillations in other channels are suppressed. The 
maximal mixing is natural in the sense that it results from the large value of the My^y^ 
matrix element relative to the diagonal M^^^^ and M^^y^ elements, without any need for fine 
tuning of the difference \My^y^ — My^y^\. 



6.4 Extraterrestrial baseline 



Finally, for very large L/E ^ ('^"^Lm/^V^) ^ km/GeV, sin^ A32 averages to | and the 
appreciable oscillations in vacuum are 



PU. 



P(iy„ 



P{v„ ^ 



sin^ 26* sin^ Aqi , 
t\ sin^ 26* sin^ Aqi , 
I - sin^ 26* sin^ Aqi 

2^3 1 
1 

2 • 



(71) 
(72) 

(73) 
(74) 
(75) 
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The solar data can then be explained if the parameters in vacuum satisfy 



SAM: 5mL„ = 5m2^ ^ m^e? , = ^ ^ ^ , (76) 



2 



LAM: 6ml^ = 6ml,c:,mh,^el + Ae',, = ^ , (77) 

VLW: 5mL„ = ^ 2m2e26ie4 , = ^ 1 , (78) 



in the three cases. 



6.5 Determination of the parameters 

In any of these scenarios in Sec. 5.1-5.3, the parameters m, €2, £3, €4, and ei are obtained 
from the data in exactly the same way, i.e., via Eqs. (|66D , ([TOD , and (|76D-(|78D. This is a 
consequence of the i^g-z/,- ambiguity. In all cases, the parameters m, €3, and £4 are related to 
the observables by 

™2 _ r 2 2 _ _ ""'■atm (^q^ 



s6i 



In the solar sector we have 



SAM: e? = ^, el = ^'^fj- , (80) 



LAM • = '^'"sun Vl-^sun g2 _ ^™Ln^3un ('g2^ 

VLW: ei~0, 62 = JfX^^ . (82) 

For the specific values ^m^^, = 2 eV^ and A^li = 2.5 x 10"^, Sml^^ = 5 x 10"^ eV^ 
and A'^fj^ = 1, and Sm'^^Y. = 4 x 10~^ eV^ and = 1 x 10"^, the model parameters are 
given m Table |. If we take instead ^m^^; = 0.3 eV^ and A'^^i = 4.0 x 10"^ (which gives the 
smallest value of Sm'j.iji allowed by the data), we get the model parameters shown in Table ^ 
In either of these two examples, the 6m'^ scale for the atmospheric neutrino oscillation can 
be adjusted simply by varying €4. Also in either case, the two heaviest masses provide relic 
neutrino targets for a mechanism that may generate the cosmic ray air showers observed 
above ^ 10^° eV |50]. We note that the model parameters in Tables H and H obey the 



hierarchies described in Eqs. (^3)"(e1 



7 Model predictions 

7.1 Resolving the Ur vs. ambiguity 

If the solar oscillations are — > z/^ as described in Sec. 5.1, then our four- neutrino model 
predicts that the atmospheric oscillations are z/^ — * u^-. On the other hand if the solar 
oscillations are z/g Vj. as in Sec. 5.2, the atmospheric oscillations are z/^ Ug. Several 
possibilities have been discussed to resolve the ambiguous assignment of u^- and as the 
oscillation partners of the z/g's in the sun and the z/^'s in the atmosphere. The Solar Neutrino 
Observatory (SNO) [^, which can measure both charge-current (CC) and neutral-current 
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(NC) interactions, will be able to test whether the solar u^s oscillate to sterile or active 
neutrinos: in the sterile case the CC/NC ratio in SNO will be unity and both CC and NC 
rates will be suppressed from the SSM predictions, while in the active case only the CC rate is 
suppressed. Of course if the CC measurement is consistent with the NC, one needs additional 
evidence to rule out the possibility that the SSM is in error. For instance, SuperKamiokande 
and SNO can also accurately measure the shape of the '^B neutrino spectrum, which would 
be distorted by oscillations. Also, a measurement of lower energy neutrinos, such as by 



the BOREXINO experiment ||52|, could also be used to detect deviations from the SSM 
spectrum. 

Turning to the atmospheric data, the possibilities to resolve the ambiguity center around 
the earth-matter effects which are possible in the I'^-i's oscillation channel but not in the 
z/^-z/,- channel; there is a relative phase difference between i/^ and Ug due to neutral current 
forward scattering, but there is no phase difference between and z/,-. The analytical analy- 
sis of matter-effects involving active and sterile neutrinos can be somewhat complicated |pO|| , 
but the Schrodinger-like evolution equations can always be solved numerically ||53|. Other 
tests have been proposed recently to resolve the z/^-z/g ambiguity. One test is to measure the 
asymmetry between downward-going and upward-coming events, for electrons and muons 
separately [Q. Various oscillation scenarios give rise to dramatically differing trajectories 
of the asymmetries versus energy for muons and electrons. The preliminary data from Su- 
perK for the individual muon and electron asymmetries suggests again that the atmospheric 
anomaly is primarily due to z/^ oscillating into either z/^ or Ug, but not z/g. By eventually 
measuring an up-down asymmetry for neutral current (NC) events (e.g. uN —>■ uNtc^), the 
ambiguity can be resolved: for the i>r case there is no NC asymmetry, whereas for the Ug case 



there is a large NC asymmetry, as shown in Ref. The ratio of the rates of NC events 
relative to the charged current (CC) events can be also used to the same end ||5^, as can 
multi-ring events ||5^. Searches for muon-less events which come from Ur, in association with 
a z/„ disappearance measurement, can also in principle distinguish between and Ug [BH]. 



7.2 New oscillation signals 

Assuming that the solar oscillations are u^. —>■ Ug, we can determine the new oscillation signals 
predicted by the model. Given the order of magnitude of the Smfj and Uaj, observable new 
phenomenology occurs for L/ E ^ 1 km/GeV in the oscillation channels 

P(z/e^z/^) ~ ^A^,^,(2-sin2A,,J , (83) 

P{ve^vr) - -^A^^l^sin^ Ktm : (84) 

where A^J^i ~ (9(1%) is the oscillation amplitude which describes the LSND results and 
^atm = l.276ml^^L/E ~ (SmltJb x 10"=^ eV^) {L/ 157km) {GeV / E) is the oscillation argu- 
ment which describes the atmospheric neutrino data. We emphasize the new predictions in 
the z/g z/^ and Vr channels: long baseline oscillations with common oscillation length 
determined by the atmospheric Aatm and common amplitude given by i times the LSND 
amplitude A'^^^. These oscillations are in addition to the z/^ z/g oscillations due to Agf,i in 
Eq. (0), which average to the value of \A'^^i in a long baseline experiment. The amplitudes 
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and lengths of these new oscillations complement the set in Table 0, which are inevitable, 
given the present data, and are therefore required in any model. 

How can the oscillation probabilities in Eqs. ( ]83| ) and (Q) be tested? A list of experi- 
ments currently underway or being planned to test neutrino oscillation hypotheses is given 
in Table ^ |^ . In each case the oscillation channel and the parameters which are expected 
to be tested are shown. 

The MINOS experiment can detect z/^ v^, or z/^ oscillations and is sensitive 

down to brn? ~ 10~^ eV^ and a mixing amplitude of 10~^, which partially overlaps the region 
of interest; see Fig. 8. If the MINOS experiment can increase its sensitivity, it will provide 
an even better test of this new phenomenology. 

Long-baseline experiments with an intense z/g or Vf> neutrino beam and which can detect 
r's can see the z/g — > oscillations in Eq. ([8^ ) and provide a definitive test of the new 
phenomenology predicted by the model. High intensity muon sources can provide simul- 
taneous high intensity z/^ and Ve (or z/^ and z^e for antimuons) beams with well-determined 
fluxes, which could then be aimed at a neutrino detector at a distant site. It is expected 
that r's will be detected through their /i decay mode and that a charge determination can be 
made, so that one can tell if the r originated from z/^ or z/g oscillations. Current 
proposals consider SOUDAN (L = 732 km) or GRAN SASSO (L = 9900 km) as the far 
site from an intense muon source at Fermilab (MC). These experiments could also observe 

— ^ oscillations via detection of "wrong-sign" muons, i.e., those with sign opposite 
to that expected from the or source. The neutrino energies are in the 10-50 GeV 
range. Assuming that low backgrounds can be achieved, the sensitivity to bw? is roughly 
proportional to the inverse square root of the detector size (given the same neutrino energy 
spectrum at the source); the brr? sensitivity does not depend on detector distance L because 
although the flux in the detector falls off with L^, the oscillation argument grows with I? for 
small 5m?L/ E. For 20 GeV muons at Fermilab and a 10 kT detector at either SOUDAN or 
GRAN SASSO, the single-event sensitivity for z/g Vt oscillations is about 8 x 10~^ eV^ 
for maximal mixing [^]. For large 5m^, the oscillation amplitude single-event sensitivity is 
roughly inversely proportional to the neutrino flux at the detector divided by the detector 
size; about 6 x 10"^ for SOUDAN and lO'^ for GRAN SASSO ig. In general, the closer 
detector has comparable bvn? sensitivity but better A sensitivity. 

The model predicts z/g oscillations with amplitude ^A^u (which ranges from 0.0006 
to 0.01) and mass-squared difference of bm^^^^ (which ranges from 3 x 10~^ to 7 x 10~^ eV^). 
The region of possible z/g — oscillations in our model and the regions which can be tested 
at the SOUDAN and GRAN SASSO sites with a neutrino beam from a high-intensity muon 
source at Fermilab are shown schematically in Fig. 8, along with the favored parameters for 
the LSND, atmospheric neutrino, and solar neutrino oscillations. Such experiments would 
be sensitive to some of the z/g — >■ region, though they may not cover the low-mass, small- 
amplitude part. These searches would also be able to test the z/g — >■ oscillations in Eq. 
(^) and the atmospheric oscillations. Additionally, long baseline experiments to 

the AMANDA detector g from Fermilab (L ~ 11700 km) or KEK (L ^ 11300 km) may 
be useful in probing oscillations with small bw? . 

If the solar oscillations are z/g u^-, then the oscillations of atmospheric neutrinos are 



J^iM and the new oscillations in Eq. (|8^) are instead z/g — Vg. Neither of these signals 



would be detectable in long-baseline experiments since the signal is z/g or z/„ disappearance 
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at the few percent level or less. The only measurable signal of the model in this case is the 
oscillations in Eq. (^). 
If the solar neutrinos oscillate into both and as given by the mixing in Eq. (p9D, any 
vacuum oscillation in Sees. 6.1-6.4 which has Vg as the final state is replaced by oscillations 
to Vs with relative probability cos^ a and to with relative probability sin^ a. Conversely, 
any oscillation which has as the final state is replaced by oscillations to with relative 
probability sin^ a and to with relative probability cos^ a. In particular, a solar z/g oscillates 
into a mixture of and with relative probability cos^ a and sin^ a, respectively, in a 
vacuum, and an atmospheric oscillates into a mixture of and with relative probability 
sin^ a and cos^ a, respectively, in a vacuum. The new oscillation signal in Eq. (|8^) for long 
baseline experiments is replaced by 

P[v, ^ V,) ~ \A^^^ cos2 a sin^ A,,^ , (85) 
P(z/e ^ V,) ~ \A^^l sin^ a sin^ A„t„ . (86) 

Also, there are new oscillations between and z/^, with the general probability for any 
baseline given by 



P{ys^Vr)= sin a cos a 2cos 6'(sin A20 + sin A30) 

+2sin2 6'(sin2 A21 + sin^ A31) - sin^ A32 - sin2 26'sin2 Aoil . (87) 



7.3 Neutrinoless double— /3 decay 

Since the M^^y^ mass matrix element is zero in Eq. (^51) , there is no neutrinoless double- 
/3 decay at tree level in our model. The present limit on M^^y^ from this process is ~0.5 
eV [Q. New experiments are under development which may measure My^y^ down as low as 
0.01 eV||H^. If a nonzero M^^y^ is found at these levels, it would be incompatible with the 
solar solutions in our models. 



7.4 Tritium decay 

If Ve is primarily associated with the lighter pair in the 2+2 model, and mi < itlq ^ m2, m^, 
then there will be no measurable effect in the endpoint of the tritium decay spectrum. Since 
only mass-squared differences are important for oscillations, the inverted 2+2 model, where 
7712,1713 -C rrii < mo, is equally viable, although it is not derivable from the mass matrix in 
Sec. 5. Then the will have an effective mass of 0.55 — 1.4 eV, which is just below the 



current limit 34 



7.5 Hot dark matter 

For m ~ 1.4 eV, approximately the largest value allowed by the LSND data, J2''t^u — 3 eV, 
which according to recent work on early universe formation of the largest structures provides 
an ideal hot dark matter component For m ~ 0.55 eV, the contribution to hot dark 

matter is much smaller. The contribution of the neutrinos to the mass density of the universe 
is given hj fly = J2 ^u/ (/i^93 eV), where h is the Hubble expansion parameter in units of 100 
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km/s/Mpc |66[; with h = 0.65 our model implies ~ 0.05. An interesting test of neutrino 
masses is the power spectrum of early galaxy sizes, to be provided by the Sloan Digital Sky 
Survey (SDSS) For two nearly degenerate massive neutrino species, sensitivity down to 
about 0.2 to 0.9 eV (depending on fl and h) is expected, providing coverage of all or part 
of the LSND allowed range (m =0.55 to 1.4 eV in our model). Also, the future MAP 



and PLANCK satellite missions, which will measure the cosmic microwave background 
radiation, should be sensitive to neutrino densities to high precision |]7D[, and in particular 
to the z/^ — > atmospheric or the z/g Vg large-angle MSW neutrino mixing solutions |[71|| . 

7.6 Resonant enhancement in matter 

The curves in Fig. 8 (in the scenario where the solar oscillations are — > z^s) assume vacuum 
oscillations in the Earth. In general, large corrections to oscillations involving z/g and Vg are 
possible due to matter. The Ve diagonal element in the effective mass-squared matrix receives 
an additional term \/2GFNf.E from its forward elastic CC interaction, and the Vg diagonal 
element receives the contribution GpNnE / \/2 (relative to the active neutrinos) because it 
does not have NC interactions. Here, A^e and Nn denote the electron and neutron number 
density, respectively. In Table ^ we display the resonant energies in the earth for the various 
vacuum 5w? values suggested by the available LSND, atmospheric, and solar data. For 
neutrinos with energies significantly above E^esi oscillations are suppressed; for neutrinos 
with energies significantly below E^esi the matter effect is negligible; at the resonant energy, 
^mat = 1 and the oscillation length is increased by 1/ V^vac- 

Some of the resonant energy values in Table ^ are of particular interest. Upcoming 
neutrinos from the atmosphere or astrophysical sources, with mass at the lower end of the 
LSND range, can have their oscillations resonantly enhanced by the earth's mantle and/or 
core. Atmospheric neutrinos below a few GeV and the SAM and LAM pp neutrinos from 
the sun also appear to be near resonance in the earth's matter. Day-night modulation of 
the solar flux due to earth-matter effects is expected to discriminate between Vj. and Vg solar 



fluxes [72, 73 1, while a precise measurement of zenith angle dependence may discriminate 



between v^- and Vg atmospheric fluxes [Q. Oscillation wavelengths commensurate with the 



size of the earth's mantle and/or core are especially sensitive. 

In our model of Sec. 5.1 with —>■ Ur atmospheric neutrino oscillations, however, these 
corrections do not significantly affect the large ml and mg mass eigenvalues as long as 
E <^ 1 TeV, and hence only modify the Sm^i oscillation argument. We have verified this 
by explicit diagonalization of the mass matrix when matter effects are included. Hence 
we conclude that the matter corrections to the mass matrix in Eq. ( ^3D probably have 
no observable consequences in all terrestrial experiments. For large L/E, such as when 
E ^ 10 MeV for solar neutrinos, the only significant effect of matter is the usual MSW 
enhancement of — »• z/^ that leads to the solar neutrino suppression of the z/g flux; in all 
other oscillation channels the matter-enhanced amplitudes are at the e^^ level or smaller. 
In the models discussed in Sees. 5.2 and 5.3, where there is a z/^- — Ug component to the 
atmospheric oscillation, matter effects as discussed here may be important in terrestrial 
experiments [EO]. 
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8 Discussion and summary 



8.1 Distinguishing the three solar solutions 

The VLW solar solution may be discriminated from the two MSW solutions by a careful mea- 
surement of the solar neutrino spectrum by SuperK and BOREXINO |75], or by determining 
the amount of seasonal variation of the ^Be and pep neutrinos , which can be measured 
by the BOREXINO experiment. The ^B neutrino spectrum as measured in SuperK and 
SNO will also be useful in discriminating between the SAM and LAM solutions The 
HERON and HELLAZ ||78| experiments would be able to measure the pp neutrino energy 
spectrum, which would also be useful in differentiating the three scenarios. 



8.2 Possible Vt—Vs mixing 

The general case with z/^-z/t- mixing is described in Sec. 5.3. The unmixed cases Eqs. ( ^5|) 
and (|5^) are obtained in the limits a — *■ and a — * ^, respectively; distinguishing between 
the unmixed scenarios is discussed in Sec 7.1. How might non-trivial mixing of Ur and Ug 
be observed, and how might the mixing angle be deduced? A mixed model would generally 
have a signature intermediate between the two unmixed signatures p^j; e.g., experiments 
measuring neutral current events for solar and atmospheric neutrinos would find a result 
between those expected for Ur and Ug. 



8.3 Summary 

An analysis of all the available data (short baseline LSND, reactor and accelerator, long 
baseline atmospheric, and extraterrestrial length solar) in terms of neutrino oscillations leads 
to the conclusion that three independent oscillation lengths are contributing. This then 
further requires mixing of at least four light-mass neutrinos. For a four light-mass neutrino 
universe, we draw the following model-independent conclusions: 

(i) the 1+3 (or 3+1) mass spectrum with a separated mass is disfavored when all the data 
(LSND, reactor, accelerator, solar, and atmospheric) are considered, leaving a spectrum with 
two nearly-degenerate pairs as preferred; 

(ii) the neutrino mixing matrix elements satisfy \Ue2/Ue3\ — \Ufj,2/U^3\ if ^m^^; ~ 0.3 eV^, 
i.e., the parameters lie near the Bugey z/g disappearance limit; 

(iii) the relation |?7^2| — l^^sl inferred from the near-maximal mixing of atmospheric z/^'s 
measured by SuperK, which together with (ii) implies \Ue2\ — jt^esl- 

Based upon the apparent need for more than three light neutrinos, we have presented 
four-neutrino models with three active neutrinos and one sterile neutrino. The models nat- 
urally have maximal z/^ — > u^- (or z/^ Ug) oscillations of atmospheric neutrinos and can 
also explain the solar neutrino and LSND results. The solar solutions can be z/g —>■ or 
z/g — z/^, and can be small- angle matter-enhanced, large- angle matter-enhanced, or vacuum 
long-wavelength oscillations; the increased statistics on the electron energy distribution and 



day/night differences of the SuperK data may further clarify the allowed regions for the 



solar solutions. The models predict z/g ^ i^r (or z/g Ug) and z/g ^ z/^ oscillations in long- 
baseline experiments with L/E ^ 1 km/GeV with amplitudes that are determined by the 
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LSND oscillation amplitude and Sm^ scale determined by the oscillation scale of atmospheric 
neutrinos. For the Uj. case, these oscillations might be seen by experiments based on 

neutrino beams from an intense muon source at Fermilab with a detector at the SOUDAN 
or GRAN SASSO sites. The i^^ — > oscillations might be seen by the MINOS experiment 
or at KEK with detectors at Kamiokande and SuperKamiokande. The models also predict 
the equality of the i/^ disappearance probability, the disappearance probability, and the 
LSND i/^ — > i/g appearance probability in short-baseline experiments. 
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Table 1: Ranges of mass- squared differences and amplitudes that provide oscillation solutions 
to the solar neutrino data within 95% CL. in the small-angle MSW, large-angle MSW and 
vacuum long-wavelenth scenarios when z/g — > Us or — > i^r The new SuperK data and 
new solar flux calculations give slightly different oscillation parameters than those quoted 
here; in particular, the (5m^ values for the VLW solution are higher, of order 4 x 10~^° eV^ 

El- 



l^e ^ 


SAM 


LAM 


VLW 


5m^ (eV^) 
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Table 2: Some typical 6m'^ and vacuum oscillation amplitudes suggested by experiment, 
and the corresponding vacuum oscillation lengths A„ = 2.5 (-E't,/GeV)((5m^/eV^)-^ km for 
representative neutrino energies. Here AU is the earth-sun distance and i?® is the earth's 
radius. 
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Table 3: Typical model parameters for the input data Smlf^i = 2 eV^, A'^^i = 2.5 x 10 
Smltm = 5 X 10-3 eV2, and = 1. 
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Table 4: Typical model parameters for the input data Smlj^i — 0.3 eV^, A^^i — 4.0 x 10 
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Table 5: Current and planned short and long baseline neutrino oscillation experiments. Stars 
denote accessible oscillation channels. The Siri^ and sin^ 29 sensitivies are given. 
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Table 6: Resonant energies E^es = 6.6 ^/T^^ {6m'^ /eV^){NjNA cm-^)"! TeV in the earth's 
core and mantle for Ve-^n/r oscillations, for some typical vacuum values for the 5m^'s and 
amplitudes suggested by the data. Here A is the oscillation amplitude and Nj^ is Avagadro's 
number. Wc have taken the core electron density to be 4.5 to 6.0 Na/ciu.^, and the mantle 
density to be 1.6 to 2.6 iV^/cm^. Resonant energies for Ug^^s oscillations are 2Ne/{2Ne — Nn) 
times larger than for ye^i^n/r oscillations, resonant energies for v^/r-i^s oscillations are 2Ne/Nn 
times larger, but resonance occurs for P rather than u since the phase difference due to matter 
has the opposite sign. 
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Figure 1: Neutrino mass spectrum, showing a possible flavor assignment for each mass 
eigenstate, and showing which mass splittings are responsible for the LSND, atmospheric, 
and solar oscillations. 
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Figure 2: Typical probability rectangles for a 2+2 model where z/2 and z/3 generate the 
atmospheric oscillations and vi and vq generate the solar oscillations. For the non-dominant 
probabilities j = 2 or 3 and A; = 1 or 0. 
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Figure 3: The three allowed two-neutrino solar solutions for u^, oscillations |45|. The 

corresponding region for z/g — >■ oscillations are similar to the z/g case. 
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Figure 4: Typical probability rectangles for some 2+2 models, (a) Rectangles deduced from 
atmospheric and CHOOZ neutrino data, with i/^ lying predominantly in the z/2 and z/3 states 
with large mixing with either v-r and/or z/^. The P^o and P^i probabilites are small, but 
the P^o and P^i probabilities are not yet determined. The partitioning of the P^o and P^i 
probabilites once the solar solution is specified are also shown for the (b) small angle MSW 
and (c) large angle MSW or vacuum long-wavelength solutions. 
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Figure 5: Typical probability rectangles for the 1+3 model when Pgs is large and P^v is small. 
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Figure 6: Constraints on P^^ and Pgs for the 1+3 model when (a) Sni^i^i = 1.7 eV^ and (b) 
Smli^i = 0.5 eV^. The LSND, Bugey, CDHS, atmospheric, and solar constraints are obtained 
by comparison of the appropriate data with Eqs. (|T^, (PPD, PT|), (P^, and (P^D, respectively. 
Not displayed is the unitarity constraint Pes + P^3 < 1. 
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Figure 7: Constraints on and (Xe defined in Eq. ( ^l]) for tlie 2+2 model wlien 5m^j; = 
1.7 eV^. Tfie LSND, Bugey, CDHS, solar and atmospheric constraints are obtained by 
comparison of the appropriate data with Eqs. (|3^) , (|29|) , (|30D , (|33|) , and (^l]), respectively. 
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Figure 8: Predicted region in the effective 5m^-sin^ 29 parameter space for z/g — *■ u-r oscilla- 
tions in the four-neutrino model (solid rectangle), which is determined by j of the LSND 



oscillation amplitude and the atmospheric neutrino z/^ 



oscillation Sm^ scale. 



z/p, oscillations from the MINOS 

oscillations 



The dotted curves show the potential limits on 
experiment ||60[ and the dashed curves show the potential limits on u, 
that can be set by neutrino beams from an intense muon source at Fermilab to detectors 
at the SOUDAN and GRAN SASSO sites for muons with energy of 20 GeV. Also shown are 
the parameters for the solar z/g — > z^s small-angle MSW oscillation. 
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